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PREFACE 


This  report  was  prepared  by  2Lt.  William  J.  Miklos  of  the  Fluid 
Mechanics  Research  Laboratory,  Aerospace  Research  .aboratories  of  Air 
Force  Systems  Command, and  Dr.  Charles  W.  Ingram,  visiting  scientist  at 
Aerospace  Research  Laboratories  employed  through  Technology  Incorporated. 
The  work  described  herein  was  requested  by  Air  Force  Armament  Laboratory 
and  was  accomplished  between  September  1974  and  February  1975.  Mr.  Paul 
D.  Shirey  of  Armament  Laboratory  served  as  program  manager. 
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SECTION  I 


INTRODUCTION 

Air  Force  Armament  Laboratory  (AFATL)  has  requested  that  the  Aerospace 
Research  Laboratories  ( ARL ) assist  in  the  aerodynamic  evaluation  of  three 
replacement  designs  for  the  MK  82  Snakeye  Bomb.  The  purpose  of  this  report 
is  to  document  a series  of  subsonic  dynamic  wind  tunnel  tests  carried  out  at 
AFATL,  transonic  dynamic  wind  tunnel  tests  conducted  at  the  Naval  Ship 
Research  and  Development  Center  (NSRDC),  and  transonic  static  tests  accom- 
plished at  Arnold  Engineering  Development  Center  (AEDC)  to  determine  the 
aerodynamic  characteristics  of  the  candidate  configurations. 

The  Snakeye  Bomb  depicted  in  Figure  1 is  an  aircraft  store  which  may 
operate  in  either  a high  or  low  drag  flight  mode.  However,  the  transonic 
static  instability  of  the  low  drag  mode  0*2,)  and  the  mechanical  unreliabil- 
ity of  the  high  drag  mode  have  necessitated  the  design  and  evaluation  of  a 
new  weapon  to  fulfull  this  dual  mode  mission  requirement.  Although  the  three 
candidate  replacements  have  this  capability,  the  mechanical  fin  decelerator 
system  of  the  Snakeye  has  been  replaced  by  a ram  air  inflatable  retarder. 

This  retarder  is  stored  in  the  stabilizer  assembly  of  the  bomb  and  is 
deployed  when  the  weapon  has  separated  from  the  aircraft. 

The  three  low  drag  candidate  configurations  which  were  tested  are 
depicted  in  Tigure  2.  While  they  consist  of  the  same  six  caliber  forebody, 
they  differ  in  the  design  of  their  tail  assemblies.  The  two  basic  configur- 
ations shown  in  this  figure  result  from  the  aircraft  carriage  requirement 
that  the  tip-to-tip  fin  span  be  no  greater  than  1.4  times  the  maximum  body 
diameter.  The  first  of  these  designs  employs  a fixed  cruciform  finned  tail 
section  similar  to  that  of  the  Low  Drag  Bomb  with  the  exception  that  it  has 
no  boattailing  and  an  increased  fin  chord.  The  other  basic  shape,  the  ISR, 
utilizes  an  inflatable  stabilizer.  The  tail  assembly  of  this  design  is  folded 
for  aircraft  carriage  and,  upon  release,  is  deployed  to  the  low  drag  config- 
uration shown  in  the  figure.  The  third  shape  simply  results  from  the  addition 
of  an  extender  to  the  second  basic  design.  All  three  configurations  will 
employ  the  same  ram  air  inflatable  decelerator  for  the  high  drag  flight  mode. 

Single  degree  of  freedom  free-osci 1 lation  wind  tunnel  tests  were 
conducted  at  AFATL  and  NSRDC  to  obtain  the  static  moment  stability  coefficient, 

Cf'iot,  and  the  damping  moment  stability  coefficient.  + Cmc"x  • These  coef- 
ficients were  determined  as  functions  of  center  of  gravity  position  for  all 
candidate  shapes  at  the  following  Mach  numbers:  0.13,  0.5  and  0.8.  The 

normal  force  coefficient  slope,  CZa,  and  the  center  of  pressure  location  of 
each  configuration  were  determineafrom  the  CMa  versus  center  of  gravity 
position  data.  Static  tests  of  the  Fixed  Fin  aesign  were  also  carried  out 
at  AEDC  to  confirm  the  results  of  the  AFATL  and  NSRDC  tests  and  to  better 
define  the  aerodynamic  characteristics  of  this  configuration  at  higher  Mach 
numbers.  Force  and  moment  measurements  at  the  flight  center  of  gravity  were 
obtained  at  Mach  numbers  of  0.3,  0.5,  0.8,  0.9,  0.95,  and  1.0. 
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SECTION  II 
Th‘EOr  Y 

During  the  AFATL  and  NSRDC  dynamic  tests,  the  model  was  supported  by  a 
vertical  strut  which  restricts  its  angilar  oscillation  to  the  yaw  plane.  As 
a result,  the  motion  was  not  affected  by  any  moments  created  if  the  center 
of  gravity  did  not  lie  below  the  strut  centerline.  The  dynamics  is  therefore 
comparable  to  a modei  freely  pitching  in  a wind  tunnel  while  supported  at 
its  center  of  gravity.  If  the  mount  contributes  negligible  friction  forces, 
the  angular  motior  is  governed  by  the  three  forces  and  their  moments  depicted 
in  Figure  3.  These  forces  are  referred  to  as  the  static  force,  damping  force, 
and  aerodynamic  lag  force  and  are  denoted  by  (5) 


(static) 

Zaa.  = Dj;  a QS 

( 1 ) 

(damping) 

■ Czq  QS 

( 2 ) 

(lag) 

Z&i  * Cz&  § CIS 

( 3 ) 

The  moments  of  these  forces  about  the  center  of  gravity 
this  case)  are  described  by  the  following  equations: 

(or  pivot  point  in 

(scatic) 

MqK  - a QSd 

‘Ot 

( 4 ) 

(damping) 

Mqq  = C‘-!q  2V  ^Sd 

( 5 ) 

(lag) 

ffca  = CM5(  QSd 

( 6 ) 

Aerobal  1 i sticians  have  resolved  that  the  damping  moment  and  the  lag  -mient 
cannot  be  measured  separately,  and  therefore  Eq s( 5 ) and  (6)  are  alw^s 
written  in  a coupled  form.  The  equations  which  are  essential  to  a single 
degree  of  freedom  stability  analysis  of  a body  with  no  aerodynamic  asym- 
metries are  therefore 

Mua  = cr-’u  “ QSd 

( 7 ) 

(Mq+Md)  u 

= (CHq*C„.)  ^ QSd 

( 8 ) 

It  is  now  possible  to  formulate  the  differential  equation  whose  solution 
will  describe  the  oscillation  of  the  body.  From  Newton's  Second  Law 


lytt  - E n ( 9 ) 


2 


Substituting  the  expressions  for  the  static  moment  and  the  damping  moment, 
one  obtains  the  desired  result, 

Iy'a  = IV  + (Mq+V  a (IQ) 


or  in  another  form 


a 


(-V')  “ 


a = 0 


(ID 


The  solution  to  this  linear,  second  order,  homogeneous  differential  equation 
is  a damped  cosine  function: 

^ a (t)  ■ a0  eAt  cos(cot+6)  (12) 


where 


and 


X='2\/'  (13) 

- J 1* 

“ ' V (14) 


Using  Eqs  (7)  and  (8)  the  above  expressions  for  the  damping 
frequency  become 


>. 


(C»-  +C>,«.) 

v "q  -a 


QSd- 

*v 


rate  and 


(15) 


and 


w = 


J C^Sd 

1 Iy 


(16) 


7 body  is  said  to  possess  static  stability  when 

C:,w  < 0 (17) 


and  is  said  to  possess  damping  stability  when 


+Cm- 


0 


(18) 
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SECTION  III 


EXPERIMENTAL  PROCEDURE 

In  both  the  AFATL  and  NSRDC  dynamic  tests,  the  free-osci llation 
technique  was  used  to  obtain  data  required  to  extract  the  desired  aerodynamic 
coefficients.  Schematics  of  the  two  test  section  support  systems  are 
presented  in  Figures  4 and  5.  The  models  were  manually  displaced  to  a high 
angle  of  attack  orientation,  and,  after  they  were  released,  the  subsequent 
motion  of  the  pointers  was  photographed  by  a high  speed  camera.  These 
photographic  data  were  then  analyzed  to  determine  the  circular  frequency 
and  damping  rate  of  the  oscillation. 

The  AFATL  subsonic  wind  tunnel  is  an  open  circuit  facility  which 
utilizes  a radial  blower  to  force  ambient  air  through  the  3x5  foot  test 
section.  While  the  maximum  attainable  velocity  is  190  feet  per  second,  all 
dynamic  tests  were  conducted  at  approximately  150  feet  per  second.  In  order 
to  insure  that  sufficient  data  were  available  to  determine  the  normal  force 
characteristics,  the  three  configurations  were  tested  with  the  pivot  point 
located  at  3.154  diameters,  3.655  diameters  (ISR  flight  center  of  gravity), 
and  3.934  diameters  (Fixed  Fin  flight  center  of  gravity)from  the  nose  of  the 
model.  The  required  transverse  moments  of  inertia  were  calculated  From 
torsional  pendulum  data  obtained  at  the  Measurements  Laboratory  of  the 
Ballistic  Aerodynamics  Research  System  of  AFATL. 

Staoility  data  at  higher  Mach  numbers  were  collected  in  the  closed 
circuit,  single  return,  5x  7 inch  wind  tunnel  of  the  NSRDC.  As  depicted  in 
Figure  5,  the  top  and  bottom  walls  are  slotted  to  allow  the  flow  to  expand 
into  a pressure  chamber  which  surrounds  the  test  section.  During  tunnel 
operation,  the  total  pressure  in  the  circuit  was  set  at  1100.0  pounds  per 
square  foot,  and  the  total  temperature  .cached  a maximum  of  98°F  at  M = u.8. 
Each  configuration  was  tested  with  the  vertical  strut  located  r. t ^x.al 
stations  which  were  3.655  diameters  and  3.934  diameters  from  the  nose.  The 
transverse  moments  of  inertia  of  the  0.625  inch  diameter  models  were  calcul- 
ated at  ARL  by  the  torsional  pendulum  method. 

Further  static  tests  of  the  Fixed  Fin  configuration  were  accomplished 
at  AEDC  in  the  four  foot  transonic  tunnel  of  the  Propulsion  Wind  Tunnel 
Facility.  Since  this  closed  circuit,  variable  density  tunnel  is  designed 
to  operate  over  a Mach  number  range  from  0.1  to  1.3,  its  test  section  is 
equipped  with  variable  pourosity  walls  to  control  flow  blockage  and  shock 
wave  interference.  Tne  one-fifth  scale  model  which  was  tested  in  the  AFATL 
tunnel  was  also  used  during  these  transonic  tests  after  it  was  modified 
to  be  compatible  with  the  existing  sting  support  equipment.  Aerodynamic 
forces  and  moments  were  measured  by  an  internally  mounted  six  component 
strain  gauge  balance  and  the  electrical  output  was  analyzed  by  an  on-line 
Ratheon  520  computer.  A more  complete  description  of  the  tunnel  and  its 
associated  equipment  may  be  found  in  Reference  6. 
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SECTION  IV 
ANALYSIS  PROCEDURE 


The  dynamic  test  data  collected  at  AFATL  and  NSRDC  were  analyzed  at  ARL 
using  a mathematical  technique  which  will  be  developed  in  this  section.  The 
method  relates  the  frequency  of  oscillation,  the  angular  damping  rate.and 
other  physical  parameters  to  the  desired  aerodynamic  coefficients.  The 
circular  frequency  of  the  motion  may  be  defined  in  terms  of  the  period  of 
oscillation. 


When  oj  is  computed  for  each  cycle,  the  values  are  averaged  and  the  result 
is  used  in  the  calculation  of  the  static  coefficient. 


A rearrangement  of  Eq  (14)  yields 

Ma  = - u)2 1 (20) 

where  the  value  of  I is  the  transverse  inertia  of  the  model  plus  the  axial 
inertia  of  the  attached  vertical  strut.  It  is  also  true  that 


■ C„aQSd 


(21) 


Equating  these  expressions  yields  the  desired  equation  which  relates  the 
angular  frequency  to  the  static  moment  stability  coefficient. 


(22) 


After  the  aerodynamic  trim  angle  and  any  pointer  misal inements  have  been 
removed  from  the  angular  data,  the  oscillation  extrema  may  be  connected  by 
a curve  which  is  described  by  the  following  equation: 


<*„,  (t)  = <>oe  (23) 

Taking  the  natural  log  of  this  expression  and  rearranging  yield 

x = in  Op,  It)  - in  uQ  (24) 

f--r0 

Therefore,  simply  calculating  the  slope  of  the  line  drawn  through  the  in  (%, 
t points  yields  the  value  of  A.  Figure  6 is  typical  of  the  plots  generated 
by  this  method.  From  Eq  (13)  it  is  found  that 


Mq  + = 21 A 


(25) 
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(26) 


Equating  the  above  expression  wit.i  Eq  (8)  yields 


X 


4 IV 
QSd2 


The  desired  coefficients  may  be  determined  from  Eqs  (22)  and  (26). 


After  Cm  has  been  determined  at  each  center  of  gravity  location,  it 
becomes  possible  to  calculate  the  value  of  £he  normal  force  coefficient 
slope  and  the  center  of  pressure  location  (''.  In  Figure  7 it  may  be  seen 
that  the  moment  with  respect  to  a pivot  at  the  axial  station  XCG  is  defined 
by  the  following  equation: 

CMQSd  = CZQS  (XCQ-XGp)  (?i\ 


or 


Cm  - cz 


XcG“XcP 

d 


(28) 


Accordingly,  the  moment  coefficient  C^*  referred  to  a pivot  center  at 
axial  station  XCG  + AX  is 

Cm*  - Cz  -XcG'TXcP  (29) 


Subtracting  Eq  (28)  from  Eq  (29)  and  rearranging  yield 


C-  = 


r *-f. 

Si Si 

AX/d 


(30) 


Since  Cz  and  CM  are  functions  of  angle  of  attack,  the  above  equation  may  be 
differentiated  with  respect  to  a: 


Cj,i*-Cm 

na  1 ‘a 

AX/d 


(31) 


where  the  moment  stability  derivatives  are  obtained  from  two  free  oscillation 
tests  with  the  pivot  located  at.  the  two  axial  stations  XlV  and  XCCi  + AX. 
Therefore,  on  a Cr,lt  versus  Xcr/d  plot,  the  slope  of  a line  drawn  through  the 
data  points  yields *C?  . 
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The  distance  between  the  pivot  point  and  the  aerodynamic  center  of 
pressure  is  now  obtained  by  using  Eq  (28)  and  differentiating  with  respect 
to  angle  of  attack: 


XCG-XCP 


Cz 


a 


(32) 


Since  the  analysis  of  the  transonic  data  was  accomplished  by  AEDC 
personnel,  a detailed  discussion  of  their  computerized  reduction  techniques 
is  neither  required  nor  desirable.  However,  it  should  be  noted  that  the 
values  of  the  static  moment  stability  coefficient  were  obtained  by  fitting 
a line  through  the  Cm  versus  a data  points  at  -2°,  0°,  and  2°  angles  of 
attack. 
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SECTION  V 


RESULTS 

Figure  8 presents  the  restoring  moment  stability  coefficient  as  a 
function  of  Mach  number  for  all  three  candidate  configurations  at  their 
nominal  flight  centers  of  gravity.  In  addition,  data  collected  by  the 
Naval  Ordnance  Laboratory  pertaining  to  the  Low  Drag  Bomb  is  plotted 
for  comparison.  Despite  the  geometrical  similarities  between  the  Fixed 
Fin  and  the  Low  Drag  Bomb,  the  Fixed  Fin  does  exhibit  less  static  stability. 
This  is  due,  in  part,  to  the  forward  movement  of  the  center  of  pressure  on 
that  design  as  a result  of  the  longer  fin  chord.  Concerning  the  Fixed  Fin 
data,  good  repeatability  has  been  obtained  between  the  NSRDC  dynamic  test 
data  and  the  results  of  the  static  tests  conducted  at  AEDC.  This  figure 
also  shows  that  the  ISR  configuration  has  a more  negative  restoring  moment 
stability  coefficient  than  the  Fixed  Fin  due  to  the  larger  fin  span  of  the 
former.  The  addition  of  the  extender  to  the  ISR  design  moved  the  fin  lift 
force  rearward  causing  a further  increase  in  the  absolute  magnitude  of  the 
static  moment  stability  coefficient. 

In  Figure  9 the  damping  moment  stability  coefficient  is  plotted  versus 
Mach  number  for  the  candidate  configurations  and  the  Low  Drag  Bomb  at  their 
flight  centers  of  gravity.  The  longer  fin  span  of  the  ISR  design  generates 
damping  moments  which  are  more  negative  than  those  of  the  Fixed  Fin  although 
this  candidate  does  possess  adequate  stability.  The  difference  between  the 
damping  coefficient  of  the  Fixed  Fin  and  the  Low  Drag  Bomb  may  be  attribut- 
able to  the  larger  fin  area  and  wider  afterbody  of  the  former.  As  may  be 
expected,  the  addition  of  the  extender  to  the  ISR  design  effects  more 
negative  pitch  damping  moment  coefficients. 

The  normal  force  coefficient  slope,  CZct,  is  presented  in  Figure  10  as  a 
function  of  Mach  number  for  the  three  candidate  configurations  and  the  Low 
Drag  Bomb.  The  slight  variation  between  the  Fixed  Fin  static  and  dynamic 
test  results  may  be  caused  by  the  different  types  of  data  acquisition  and 
analysis  techniques  employed.  The  less  exact  method  used  to  extract  the 
coefficient  from  the  dynamic  data  is  believed  to  cause  a significant  portion 
of  the  difference  exhibited.  As  shown  in  the  figure,  the  ISR  configurations 
have  more  negative  normal  force  coefficient  slopes  than  the  Fixed  Fin  design. 

The  center  of  pressure  locations  of  the  four  bomb  configurations  as  a 
function  of  Mach  number  appear  in  Figure  11.  The  effect  of  the  longer  fin 
chord  of  the  Fixed  Fin  design  is  again  evident  as  this  configuration  exhibits 
the  most  forward  center  of  pressure.  Because  of  the  addition  of  the  ex- 
tender, the  center  of  pressure  of  the  ISR-E  is  more  rearward  than  that  of 
the  ISR. 
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SECTION  VI 


CONCLUSIONS  AND  RECOMMENDATIONS 

Static  and  dynamic  wind  tunnel  tests  of  the  Fixed  Fin,  ISR,  and  ISR-E 
candidate  Snakeye  Bomb  replacements  have  shown  that  the  normal  force  coef- 
ficient slope  and  the  static  and  damping  moment  stability  coefficients  of 
the  ISR  designs  are  more  negative  than  those  of  the  Fixed  Fin.  However, 
this  variation  in  magnitude  of  the  aerodynamic  stability  coefficients  is  not  a 
sufficient  criterion  for  eliminating  any  design  as  a potential  replacement. 
Additional  Magnus  and  roll  stability  tests  should  be  accomplished  and  the 
results  utilized  in  analytical  stability  evaluations  before  a decision  is 
made. 


Figure  1.  MK  82  Snakeye  Bomb 
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Figure  4.  AFATL  Model  Support 
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Figure  5.  NSRDC  Transonic  Test 
Section 
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Figure  7.  Nomenclature  for  C7^  Determination 
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LIST  OF  SYMBOLS 


Pitching  moment  coefficient 

JL 

QSd 


-M 


Static  moment  stability  coefficient  (radian  l) 

= !£«  =VL 

a 3a  aQSd 


CM 


Aerodynamic  lag  moment  stability  coefficient  (radian  J) 

r . 3C„  . M44 

'"Md  ~ ,4d  ■ 44,., 

3oT7  ottQSg 


2V 


2v 


Damping  moment  stability  coefficient  (radian-1) 
r = = 

Mq  HflSd 


Normal  Force  Coefficient 


Normal  Force  Coefficient  slope  (radian  *) 

r,  = _ ^aa 

'a  3a  aQS 


Diameter  of  bomb  forebody  (feet) 

Total  moment  of  inertia  (slug  foot2) 

* " *y  + *x  mount 

Transverse  moment  of  inertia  of  model  (slug  foot  ) 
Pitch  moment  (foot  pounds) 


M = CMQSd 

Restoring  moment  stability  derivative  (foot  pound/radian) 

Ma  - Cm  QSd 

a 

Aerodynamic  lag  moment  stability  derivative  (foot  pound 
second/ radian) 

"a  ■ V jQSd 
Y\i 
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LIST  OF  SYMBOLS  (continued) 


Damping  moment  stability  derivative  (foot  pound  second/radian) 
Mq  = cMq  2$QSd 

Pitching  angular  velocity  (radian/second) 

Dynamic  Pressure  (pound/ foot2) 

Q - i P V2 


Frontal  area  (feet2) 

S = nd2/4 
Time  (seconds) 

Period  of  oscillation  (seconds) 

Total  velocity  (feet/ second) 

Distance  measured  along  centerline  of  body  (feet) 

Static  force  per  radian  of  angle  of  attack  (pound/radian) 

Aerodynamic  lag  normal  force  per  radian/second  of  angular  rate 
(pound  second/radian) 

ZA  - 3QS 

Damping  normal  force  per  radian/second  of  angular  rate 
(pound  second/radian) 

zq  = czq  2^s 


Greek 


Angle  of  attack  (radians) 

Maximum  angle  of  attack  (radians) 
Initial  angle  of  attack  (radians) 
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LIST  OF  SYMBOLS  (continued) 


Phase  angle  (radians) 

Damping  rate  of  single  degree  of  freedom  motion  (second-1) 

Air  density  (slug/ foot3) 

Angular  rate  of  single  degree  of  freedom  motion  (radian/second) 

Derivatives 

First  time  derivative  (second-1) 

<•>■&:(> 

derivative  (second-2) 

) 


Second  time 
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